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Abstract 

If the X(3872) is a loosely-bound D*^D^ /D^D*^ molecule, its inclusive production rate can be 
described by the NRQCD factorization formalism that applies to inclusive quarkonium production. 
We argue that if the molecule has quantum numbers J^'-^ = l"*"*", the most important term in 
the factorization formula should be the color-octet ^Si term. This is also one of the two most 
important terms in the factorization formulas for Xcj- Since the color-octet ^Si term dominates 
Xcj production for many processes, the ratio of the inclusive direct production rates for X and Xcj 
should be roughly the same for these processes. The assumption that the ratio of the production 
rates for X and Xcj is the same for all processes is used to estimate the inclusive production rate 
of X in i? meson decays, decays, and in pp collisions. 

PACS numbers: 12.38.-t, 12.38.Bx, 13.20.Gd, 14.40.Gx 
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The X(3872) is a mysterious charmonium-like state discovered by the Belle collaboration 
in 2003 pj. It was discovered through the exclusive decay XK^ followed by the 

decay X ^ J/iIj tt+tt". The discovery was confirmed by the CDF collaboration which 
observed the inclusive production of the X meson in pp collisions. The X has also been 
observed in the discovery mode by the BaBar collaboration 3] and in pp collisions by the 
DO collaboration Q. 

The interpretation of the X meson remains a mystery 0, 0, H, S HH EM E^; E^ ■ 

The decay mode J/?/' 7r"'"7r~ suggests that it is a charmonium state, but its decay prop erties 
seem to be incompatible with its identification with any charmonium state jl^ . [l4 llq l . The 
Belle collaboration found that the width of X is small compared to that of other charmonium 
states above the DD threshold: Vx < 2 MeV at the 90% confidence level [ij. In addition 
to the discovery decay mode J/ip vr+vr", the X(3872) has been observed in the decay modes 
J/ip 7r+7r~7r° and J/ip 7 j^. Upper limits have been placed on the branching fractions 
for several other decay modes of the X, including D^D^, D^D^, D^D^ti^ [l^, Xcil, Xc27, 
J/^n^ir^ and J/ip7] 0. 

Another possibility is that the X meson is a D*°D^/D^D*^ molecule BHElEl- This 
possibility is motivated by the fact that its mass is extremely close to the D*^D^ threshold. 
The mass obtained by combining the measurements of the four experiments is Mx = 3871. 9± 
0.5 MeV 0. Subtracting the D*^D^ threshold energy, we obtain Mx - Moo - Md*o = 
+0.6 ± 1.1 MeV. For this molecule to have a substantial branching fraction into J/ip vr+vr", 
it must be a bound state, which requires that its binding energy Ex = Mj^o + Md*o — Mx be 
positive. (Note that the central value of Ex from the mass measurement is negative.) If the 
binding of the D*^D^ molecule arises from pion exchange, the most likely channel is S-wave 
with quantum numbers J^'^ = . The observed binding energy Ex is small compared 
to the natural energy scale associated with pion-exchange interactions: Ex -C m^/2yU ~ 10 
MeV, where fi is the reduced mass of the Z)*° and D^. The small binding energy implies 
that the D*^D^ scattering length a is large compared to the natural length scale l/m,r ~ 1-4 
fm. As a consequence of the large scattering length, the molecule has universal pro perties 
that are determined by a and insensitive to the shorter length scales of QCD [ll|. The 
simplest example is the binding energy itself: Ex = l/{2^a?). Another exarnple is the 
rate and momentum spectrum for the decays X —>■ D^D^tt^ and X — > D^D^'j [9]. A third 
example is the rate for the decay T(45') —>■ Xh^h^, where h"^ and h" are light charged 
hadrons (2o|. Another universal feature associated with the large scattering length is an 
enhancement in the D*^D^ invariant mass distribution near the D*°D° threshold in the 
decay B+ D*^D^K+ 

A challenge to any interpretation of the X(3872) is to understand its production rate. 
The Belle amd Babar collaborations have measured the product of the branching fractions 
for the two decays associated with its discovery 

Br[5+ ^ Xir+] Br[X ^ J/ip tt+tt"] = (1.3 ± 0.3) x 10"^ (1) 

The only other information about the production of X that is available is that its inclusive 
production has been observed in pp collisions. 

In this paper, we will discuss the inclusive production of the X(3872). We will point out 
that if the X meson is a loosely-bound S-wave 0*^0^/0^0**^ molecule, its inclusive pro- 
duction rate can be analyzed using the NRQCD factorization formalism [22] . We will argue 
that if its quantum numbers are 1+"*", the most important term in the NRQCD factorization 
formula is the color-octet ^5*1 term. This is also one of the two most important terms in the 
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NRQCD factorization formulas for XcJ- Thus for all processes in which XcJ production is 
dominated by the color-octet '^Si term, the ratio of the inclusive production rates for X and 
XcJ should be roughly the same. The assumption that the ratio of the inclusive production 
rates for X and Xcj should be the same will be used to estimate the inclusive production 
rate of X in B meson decays, decays, and in pp collisions. 

The production of the X meson involves many length scales. The creation of the cc 
pair with small relative momentum requires a hard-scattering process with a short- distance 
length scale l/rric- The evolution of the cc into a color-singlet hadronic system, either 
a charmonium or a pair of charmed mesons, occurs over a longer length scale l/rricV or 
l/rrtcv'^. The evolution of the charmed mesons into one of the four possibilities DD, D*D, 
DD*, and D*D* occurs over an even longer length scale Finally the binding of D*^D^ 

into the molecular state X occurs over the very long length scale a. We will use the NRQCD 
factorization formalism to exploit the separation between the short-distance scale rric and 
all the longer distance scales of order l/rrirV and larger. 

The NRQCD factorization formalism 22| is used to describe the inclusive production of 



charmonium. This formalism exploits the separation of scales between the hard-scattering 
process that creates the cc pair and the hadronization process by which the cc pair evolves 
into a color-singlet charmonium. The NRQCD factorization formula expresses the inclusive 
cross section for producing the charmonium state if as a sum of products of cc cross sections 
and NRQCD matrix elements: 

a[H + anything] = ^ a[cCn + anything] (O^). (2) 

n 

The cc cross section a is the inclusive cross section for producing a cc pair in the color 
and angular-momentum channel labeled by n. It involves momentum scales of order rric and 
larger. The corresponding NRQCD matrix element (O^) is the vacuum expectation value of 
an operator that creates a cc pair at a point in the channel specified by n, projects onto states 
that in the asymptotic future include the charmonium H, and then annihilates the cc pair at 
the creation point in the channel specified by n. The NRQCD matrix elements involve only 
momentum scales of order rricV or smaller, where v is the typical relative velocity of the charm 
quark in the charmonium. The sum in Q includes infinitely many terms, so it is useful only 
if the sum over channels can be truncated. The truncation can be justified by the velocity- 
scaling rules for quarkonium, which specify how each of the NRQCD matrix elements scale 
with V. For P-wave charmonium states, there are two NRQCD matrix elements at leading 
order in v: a color-singlet P-wave matrix element and a color-octet S-wave matrix element 
j23| . If the NRQCD expansion is truncated at this order, the inclusive cross section for Xcj 
can be written 

a [xc J + anything] ~ afcci (^Pj) + anything] {0^''{^Pj)) 

+ a[ccsCSi) + anything] {O^^^ {''Si)). (3) 

The subscript on cc and on the operator O indicates the color channel and the argument in 
parentheses indicates the angular-momentum channel. The relative sizes of the two terms 
in Eq. (jS)) depends on the NRQCD renormalization scale [i^. A change in the NRQCD 
renormalization scale for the operator Oi{^Pj) can be compensated by adding the operator 
Osi^Si) with a coefficient of order a^. For processes in which a color-octet cc pair can 
be created at leading order in a^, the color-octet ^5*1 term in the cross section is likely to 
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dominate for a suitable choice of the renormahzation scale. Heavy quark spin symmetry 
implies that the matrix elements in Eq. Q are proportional to 2J + 1 [22]. Thus the 
factorization formula can be written 

a[xcj + anything] ~ (2 J + l)a[cci{^Pj) + anything] (Of "(^Pq)) 

+ (2J + l)a[cc8('5i) + anything] {Of\^Si)). (4) 

The only dependence on J in the coefficient of (Of '(^^i)) is the explicit factor of 2 J + 1. 
Thus processes that are dominated by the color-octet '^5'i term can be identified experimen- 
tally by the cross sections for Xci and Xc2 being in the proportions 3:5. 

The NRQCD factorization formalism applies equally well to a pair of charm mesons, such 
as D*^ and Z)°, that are produced with a relative momentum of order ttIcV or smaller. For 
example, the inclusive cross section for producing a D*^D^ pair with kinetic energy in its 
rest frame less than -Emax? where -Emax ^ "^c? can be written 

a[D*5(E^ax) + anything] = ^l^^n + anything] 

n 

We have used D*D{Em^^) as a short-hand for D*^D^ states with kinetic energy less than 
-f'max- The cc cross sections a are the same as in In the matrix elements, the NRQCD 
operators include a projection onto states that in the asymptotic future include a D*^D^ 
pair with kinetic energy less than -Emax- If the X meson is a loosely-bound D*^D^ / D^D*^ 
molecule, the NRQCD factorization formula can also be applied to its inclusive production 
rate. The inclusive cross section for producing X can be written 

cr[X -|- anything] = a [cCn + anything] {O^). (6) 

n 

In the matrix elements, the NRQCD operators include a projection onto states that in the 
asymptotic future include an X(3872). 

The NRQCD factorization formulas in Eqs. (0) and © are useful only if the expansions 
can be truncated. Unlike the case of quarkonium, there are no velocity scaling rules governing 
the NRQCD matrix elements for D*^D^ or X that can be used to justify the truncations. 
The velocity-scaling rules would be applicable if the cc pair evolved into a heavy hadron H 
with constituents c and c only through the emission of real or virtual gluons from the c or 
c. Velocity-scaling rules do not apply if the evolution involves the emission of real or virtual 
gluons from light quarks or antiquarks in the hadron H. Thus the presence of u and u as 
valence quarks in the D*^ and guarantees that the velocity scaling rules are inapplicable. 

Another class of cc systems for which the velocity scaling rules are inapplicable is charmo- 
nium hybrid states, whose valence content is ccg [2^. Chiladze, Falk and Petrov pointed out 
that the NRQCD factorization formula could nevertheless be applied to the inclusive pro- 
duction of charmonium hybrid states. We can use similar arguments to justify the NRQCD 
factorization formula for the inclusive production of X. The NRQCD matrix elements in 
the NRQCD factorization formula for inclusive X production in Eq. (jHI) are proportional 
to the probabilities for the evolution of the cc pair into a color-singlet ccuu system, its 
subsequent evolution into an S-wave D*^D^ or D^D*^ system, and finally the binding of 
the charm mesons into the molecular state X. Since the X is weakly bound, the charm 
mesons are almost at rest. The momentum scale of the c and c in the X is therefore that 
of the c and c in the charm mesons, which is of order Aqcd- This is small compared to the 
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momentum scale of the hard-scattering process that creates the cc pair, which is of order rric 
or larger. This guarantees that the hard scattering amplitude can be expanded in powers 
of the relative momentum of the cc. The convergence of this expansion must come from a 
suppression of the NRQCD matrix elements of higher dimension operators. The convergence 
suggests a hierarchy of NRQCD matrix elements according to the dimension of the NRQCD 
operators. In particular, matrix elements associated with cc channels with higher orbital 
angular momentum quantum numbers should be suppressed. This motivates a truncation to 
a few of the lowest dimension operators in each of the four color/spin channels: color-singlet 
spin-singlet, color-singlet spin-triplet, color-octet spin-singlet, and color-triplet spin-triplet. 
The most extreme truncation of the NRQCD factorization formula would be to include only 
the four S-wave matrix elements: (C»f (^So)), {0^{^Si)), {0^{^So)), and (Cf (^^i)). 

Voloshin has pointed out that if the X(3872) is a D*^D^ molecule with J^^ = the cc 
pair in the X is necessarily in a spin-triplet state jiij. Heavy-quark spin symmetry implies 
that X(3872) is produced predominantly from cc pairs that are created in spin-triplet states. 
Thus the spin-singlet matrix elements for the production of X(3872) are suppressed. If the 
NRQCD factorization formula is truncated to include only the S-wave matrix elements, the 
leading terms are 

(t[X + anything] ^ a-[cci(^Si) + anything] (Of(^^i)) 

+ a[ccsCS,) + anything] {OiCSi)). (7) 

We now argue that the color-octet term in Eq. ((7j) will often dominate. First we note that 
since the constituent charm mesons in the X are color-singlets, the probabilities for the cc 
pair in the X to be in a color-singlet and color-octet state are 1/9 and 8/9, respectively. The 
naive expectation is that the probabilities for forming an X from the color-singlet cc pair 
and from each of the 8 color-octet cc pairs should be equal. Thus unless the color-singlet 
parton cross section a is larger than the color-octet parton cross section by a factor large 
enough to compensate for the color factor of 8, we would expect the color-octet term in 
Eq. (j?!) to dominate. This reduces the NRQCD factorization formula to a single term: 

a[X + anything] ^ a[ccsCS,) + anything] {O^CS,)). (8) 

We have argued that the color-octet ^5*1 term should be the most important term in the 
NRQCD factorization formula for the inclusive production of X, in which case it reduces 
to Eq. (jHJ. The color-octet ^5*1 term is also one of the two leading terms in the NRQCD 
factorization formula for Xcj, which is given in Eq. (jS)). For many production processes, the 
color-octet ^Si term in the NRQCD factorization formula for Xcj will dominate for a suitable 
choice of the NRQCD renormalization scale. For these processes, the ratio of the inclusive 
cross sections for X and Xcj will be the ratio of the color-octet ^Si matrix elements: 

a[X + anything] ^ (OjeS,)) 

a[xcj + anything] ^ (2 J + l){Of\^Sr)) ' ^ ^ 

The cross sections in Eq. should be interpreted as those for direct production, which 
means that the feeddown contributions from decays of heavier cc states, such as ip{2S) — *■ 
XcJl, have been subtracted. We will use the assumption that Eq. Q holds for most pro- 
duction processes to estimate the inclusive production rates of X in various processes. One 
could ignore the motivation based on NRQCD and simply take Eq. © as a phenomenological 
assumption. 
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We first consider the inclusive production of X(3872) in B meson decays. Tliis decay 
produces a mixture of B~^, and their antiparticles. The inclusive branching fraction 
for the direct production of xa from the decays of this mixture of B mesons has been 
measured to be (3.3 ± 0.5) x 10^^. (The adjective "direct" means that the feeddown from 
the decay ipi^S) Xcil has been removed.) The inclusive branching fraction is the sum 
of the branching fractions for the 2-body decay modes XciK and XciK*, the 3-body decay 
modes, etc. The fraction of the inclusive decay rate that comes from the 2-body decay mode 
XciK is (16 ± 6)%. The fraction of the inclusive decay rate that comes from 2-body decays 
is determined largely by the available hadronic energy, which does not vary dramatically 
among the charmonium states. For example, the 2-body decay B —>■ J/ipK accounts for 
(11.8 ± 1.5)% of the direct inclusive rate for B —>■ J/tp + anything. The 2-body decay 
B — > %1){2S)K accounts for (18 ±4)% of the inclusive rate for B — > ip{2S) -f- anything. These 
fractions all differ by less than a standard deviation. We will use the central value of the 
observed fraction for Xci as an estimate of the fraction for X: 



Br[i? — > X + anything] B^^^[B Xci + anything] 

Combining this with the Belle measurement in Eq. (^, we obtain the estimate 

By[B ^ X + anything] Br[X J/^ip tt+tt^] ^^x 10"^ (11) 

After multiplying by the 12% branching fraction of J/ip into leptons, the product of the 
three branching fractions is about 10~^. With the Belle and Babar collaborations having 
each accumulated more than 10^ BB pairs, the predicted rate is large enough that the 
inclusive production of the X should be observable in the data from the B factories. 

Dividing the estimate in (jllj) by the measured inclusive branching fraction for B decay 
into Xcii we obtain the product of the constant in (jU)) and the branching fraction for the 
decay X ^ J/ip tt+tt^: 

a[X + anything] ^^^^ ^ ^ ^ ^^2) 



c^dirlXci + anything] 

This can be used to estimate the inclusive production rate of X in any process for which 
the inclusive production rate of Xci has been measured. 

In the decay of the Z^, charmonium is produced primarily by the decay of the into bb, 
followed by the fragmentation of the 6 or 6 into a heavy hadron, and then the decay of the 
heavy hadron into charmonium. The inclusive branching fraction of Z^ into a charmonium 
if is to a good approximation the sum over b hadrons of the inclusive branching fraction for 
the b hadron to decay into H, weighted by the probability for a b quark to fragment into the b 
hadron, and multiplied by 2 to account for the b and the b. The inclusive branching fraction 
for the decay of Z^ into Xci has been measured. The direct branching fraction obtained by 
subtracting the small feeddown from ip{2S) decays is (2.7 ± 0.8) x 10~^. Multiplying this 
branching fraction by the number in (|12p. we obtain the estimate 

Br[Z° ^ X + anything] Br[X ^ J /i> tt+tt-] ^ 7 x 10"^ (13) 

After multiplying by the 12% branching fraction of J/ip into leptons, the product of the 
three branching fractions is about 10~^. With each of the four detectors at LEP having 
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collected millions of hadronic Z'^ decays, the predicted rate is large enough that the X may 
be observable in the LEP data. 

We next consider the inclusive production of X(3872) in pp collisions at the Tevatron. 
The inclusive sample of X mesons at the Tevatron consists of two components: a 6-decay 
component that comes from inclusive decays of b hadrons into X, and a prompt component 
consisting of X mesons that are produced directly in the pp collisions. The two components 
can be separated by using a vertex detector. Those X mesons for which the tracks of the 
pions from the decay X —>■ J/ijj 7r~^n~ and the leptons from the decay J/ip ^ V^i' can be 
traced back to the pp collision point are mostly prompt, while those for which the tracks of 
the charged particles can be traced back to a secondary vertex come from h decay. 

In Run I of the Tevatron, which consisted of pp collisions at center-of-mass energy of 1.8 
TeV, the CDF collaboration measured the fraction of prompt J/ip^s that come from Xc decay 
and the ratio of the prompt cross sections for producing Xci oi Xc2- The fraction of prompt 
J/ip^s with transverse momentum px > 4: GeV and pseudorapidity satisf ying \rj\ < 0.6 that 
come from the decays of Xci or Xc2 is 0.297 ± 0.017(stat) ± 0.057(syst) j26|. The ratio of 
the prompt cross sections for producing Xci and Xc2, where the Xc subsequently decays into 
J/ip 7 with the J/^j having < 4 GeV, is 1.04 ± 0.29(stat) ± 0.12(syst) HJ. Combining 
these results with the estimate in (fT^ . we obtain the estimate 

a[X + anything] ^^^^ _^ ^ ^ ^ ^^_3^ ^^^^ 



crlJ/ip + anything] 



The left side of Eq. fll4j) is simply the fraction of prompt J/ip^s that come from decays 
of X multiplied by the ratio of the branching fractions for X ^ J/ip 7r"'"7r~ and X 
J/ip +anything. The fraction of prompt J/ip^s that come from decays of ■0(25') are observed 
to increase with pr from (7 ± 2)% at pt = 5 GeV to (15 ± 5)% at pr = 18 GeV 0. In 
contrast, the fraction of J/ip^s that come from decays of Xc seemed to decrease shghtly with 
Pt ii^ • In Run II of the Tevatron, improvements in the detectors will allow these fractions 
to be measured as a function of pr down to pt = 0. Our assumption (0) implies that the 
fraction of prompt J/ip^s that come from decays of X should have the same pT-dependence 
as the fraction of J/ip^s that come from decays of Xcj- 

We have pointed out that the NRQCD factorization formalism that describes inclusive 
charmonium production can also be applied to the inclusive production of pairs of charmed 
mesons with sufficiently small invariant mass. If the X(3872) is a D*^D^ /D^D*^ molecule, 
the NRQCD factorization formalism also applies to the inclusive production rate of X. We 
argued that if X has quantum numbers 1^^, the leading term in the NRQCD factorization 
formula for X should be the color-octet ^5*1 term. We used the assumption that the ratio 
of the inclusive production rates of X and Xcj should be roughly the same for all processes 
to estimate the inclusive production rate of X(3872) in various high energy processes. Mea- 
surements of inclusive production rates consistent with these estimates would be consistent 
with the interpretation of the X(3872) as a D*°D^/D^D*^ molecule with J^'^ = 1++. 

I thank G.T. Bodwin for a useful discussion. This research was supported in part by the 
Department of Energy under grant DE-FG02-91-ER4069. 
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